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Abstract 

We calculate the spectra of produced thermal photons in Au+Au collisions taking into account 
the nonequilibrium contribution to photon production due to finite shear viscosity. The evolution 
of the fireball is modeled by second-order as well as by divergence- type 2+1 dissipative hydrody- 
namics, both with an ideal equation of state and with one based on Lattice QCD that includes 
an analytical crossover. The spectrum calculated in the divergence-type theory is considerably 
enhanced with respect to the one calculated in the second-order theory, the difference being en- 
tirely due to differences in the viscous corrections to photon production. Our results show that the 
differences in hydrodynamic formalisms are an important source of uncertainty in the extraction 
of the value of n/s from measured photon spectra. The uncertainty in the value of n/s associated 
with different hydrodynamic models used to compute thermal photon spectra is larger than the 
one occurring in matching hadron elliptic flow to RHIC data. 
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I. INTRODUCTION 



Experiments of ultrarelativistic heavy ion collisions performed at BNLs Relativistic Heavy 
Ion Collider (RHIC) create a hot and dense medium of quarks, antiquarks and gluons called 
the quark-gluon plasma (QGP) [l-4]. The understanding of the transport properties of 
the QGP, and the nature of confinement and the phase transition from this system to 
interacting hadrons is a central topic of modern high-energy physics. One of the most 
important discoveries at RHIC is the large elliptic flow in non-central Au+Au collisions, 
which is a clear indication of collective behavior. This implies that the QGP has a very low 
viscosity-to-entropy ratio rj/s not much larger than the AdS/CFT lower bound 1/4tt Jsj]. 
By now, it is generally agreed that the QGP thermalizes on times < 2.5 fm/c and that it 
behaves as a fluid with one of the lowest viscosity-to-entropy ratios ever observed in nature 
r]/s < 0.5 [l jj 6 8]. 

The evolution of the fireball created at RHIC has been described efficiently using relativis- 
tic hydrodynamics [l-4, 9]. Ideal hydrodynamics has been partly successful in exp laining the 



observed collective flow at low transverse momentum and in central collisions [10| . However, 
it is difficult to fit the data with it when a realistic equation of state (EoS) is used; see [3|,|4]. 
The relativistic generalization of the Navier-Stokes equation is plagued with causality and 
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stability problems |9j, so one must use the so-called second-order theories (SOTs) 
These theories rely on an expansion of the viscous tensor in velocity gradients, neglecting all 
orders higher than the second. They are unreliable in situations where these gradients are 
strong, and indeed they are known to fail, for example, in the description of strong shocks 
[ljj. It is then valuable to develop alternative theories, not limited to weak velocity gra- 
dients, to provide at least an estimate of the expected accuracy of the gradient expansion. 
With these aims in mind, in [ijj one of the authors and E. Calzetta developed an hydrody- 
namical description of a conformal field within the framework of the divergence-type theories 



(DTT) developed by Geroch [15J (see also [16|, ll7|]). DTTs do not rely on velocity gradient 
expansions and therefore go beyond second-order theories. This formalism was then applied 
in [8| to study Au+Au collisions, obtaining results in good agreement with SOTs and with 
data on elliptic flow. It was found that the nonequilibrium correction to the particle dis- 
tribution function (which is obtained from Grad's quadratic ansatz) is considerably larger 
in the DTT. This fact introduces significant uncertainty in the values of rj/s that can be 
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inferred by matching the result of different dissipative hydrodynamic theories to hadron 
multiplicity, < pt > and elliptic flow data. 

In this paper we focus on the spectra of thermal photons which are produced during the 
evolution of the fireball created in Au+Au collisions at ^s^n = 200 GeV. Thermal photons 
are produced during the entire space-time evolution of fireball, but, since they participate 
only in electromagnetic interactions, they have a larger mean free path compared to the 
transverse size of the hot and dense matter created in nuclear collisions. Therefore, photons 
created in the interior of the plasma pass through the surrounding matter without any 
interaction, providing information on properties of bulk matter and not only on its surface. 
For these reasons, the emission of photons has become a useful signature of the QGP and is 



currently being intensively studied. See 
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Thermal photon spectra have been studied within the framework of ideal hydrodynamics 

2J25I. 



18| and references therein) and using the Israel- Stewart formalism 



Recently, 



Baeuchle and Bleicher 26] have calculated photon spectra using a hydro-kinetic hybrid ap- 
proach in combination with the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) 
transport model in 3+1 dimensions j^. Moreover, the radial and elliptic flows, as well as 
the spectra and the correlations of photons are becoming a useful tool to measure or to 
further constrain the value of ri /so 



by several recent works 
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the nuclear matter created at RHIC, as has been shown 

Our purpose here is to compare the thermal photon spectra obtained using the SOT 
and the DTT 2+1 hydrodynamics to model the evolution of the fireball, both for an ideal 
and a realistic EoS which includes a QGP-hadron analytic crossover as suggested by Lattice 
QCD calculations. To this aim, we calculate the spectra of produced thermal photons 
in Au+Au collisions considering the processes of Compton scattering, gg-annihilation and 
bremsstrahlung in the QGP phase, and tttt — > pry, np — > 7T7 and p — > tttcj in the hadron 
phase. 

The nonequilibrium contribution to photon production due to finite shear viscosity (both 
in the QGP and hadronic phases) is taken into account through Grad's quadratic ansatz. 



Very recently Bhatt, Mishra and Sreekanth 



25] have calculated photon spectra including 



viscous corrections during the QGP phase using Israel-Stewart hydrodynamics and taking 



into account both shear and bulk viscosity, while Dusling 



2l| has calculated the nonequi- 



librium correction to photon production due to Compton scattering and qq annihilation at 
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leading-log order. We will compare our results to those of |2l|, |25| later on. 

We find that the use of a realistic equation of state significantly enhances the photon 
spectra, in line with the results of previous calculations (see for instance 18]). The spectrum 
calculated in the DTT turns out to be considerably larger than the one calculated in the SOT. 
The difference in the spectra calculated in both formalisms is entirely due to differences in the 
nonequilibrium correction to photon production, which is considerably larger in the DTT. 
Our results indicate that differences in hydrodynamic formalisms constitute an important 
source of uncertainty in the precise determination of r]/s from observables, such as photon 
spectra, that depend on the nonequilibrium distribution function. 

It has been shown before 0, [y] that the nonequilibrium correction to the distribution 
function has also a strong influence on hadronic observables. The calculation of photon 
and hadron observables using different hydrodynamic models and its matching to data may 
therefore provide a way of constraining the form of the nonequilibrium correction to the 



Ml- 



distribution function and thus help improve the description of the freeze-out process 

We note that in our calculations we neglect bulk viscosity, which is known from Lattice 
QCD simulations to become significant near the critical temperature [32]. It has been 



shown in 



33] that bulk viscosity influences the space-time evolution of the fireball created 



in heavy-ion collisions, thus modifying the thermal photon spectra and increasing photon 



production 



23]. The distribution function gets an additional correction coming from the 



bulk viscosity, which also modifies the spectra and puts severe limitations to the reliability 



of Grad's quadratic ansatz 25 



31] . Moreover, we do not take into account prompt photons 



from hard scatterings of partons in the initial nuclei 18|, |34j nor jet-medium photons 29] 
and we only consider thermal photons. For these reasons, we do not attempt to carry out a 
comparison with RHIC data [20] , but we focus instead on the comparison between different 
viscous hydrodynamic formalisms and on the influence on thermal photon spectra of the 
viscous correction to the distribution function. 

The paper is organized as follows. In section [III we briefly describe the hydrodynamic 
formalisms used to model the fireball evolution, and then describe the processes taken into 
account in the calculation of photon spectra, including viscous corrections. We present and 
discuss our results in section 1111} and conclude in section [TV] 
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II. THEORETICAL SETUP 



A. Evolution of the fireball 

In this section we present the hydro dynamic equations of the SOT and of the DTT for 
boost-invariant flow in 2 + 1 dimensions, for a conformal fluid. We will consider the SOT 



developed 12J, which is based on conformal invariance and generalizes the Israel- Stewart 



formalism. We employ Milne coordinates defined by proper time r = \/t 2 — z 2 and rapidity 
ij) = arctanh(2;/t). The fluid velocity is u = (u T , u x , u v , 0) and is normalized as = 1. 
The stress-energy tensor for dissipative relativistic hydrodynamics is 

T" v = puV - pA"" + W u with 

(1) 

A"" = - v?u v 

where p and p are the energy density and the pressure in the local rest frame, and IP^ is the 
viscous shear tensor which is transverse (u fl U tlu = 0), traceless and symmetric. The tensor 
A^ is the spatial projector orthogonal to For a conformal fluid we have Tjf = 0, so 
p = 3p and the bulk viscosity vanishes. 

In what follows, Latin indices stand for transverse coordinates (x,y), is the geomet- 
ric covariant derivative, D = u^D^ and V M = A^D^ are the comoving time and space 
derivatives, respectively, are the Christoffel symbols and < . . . > denote the spatial, 
symmetric and traceless projection of a tensor. The hydrodynamic equations read 

(p + p)Du i = c 2 s (g^d lP - u l u a d a p) - A^IT* 

^ = -( P +p)v/ + n% 



where c 2 s = dp/ dp is the speed of sound, = V <fl u u> is the first-order shear tensor, and 

(3) 



7/ ■ ?/ ■ 

DgU aP = TT a d T — + — d T U ia + dtlT 



In the SOT, the evolution of the shear tensor is given by 



12| 



d T Ii ia = — — IT a VX —IF* + -^—a ia 

3u T r n u T r 7r u T 

^— U< l U a> ^ £ ^Du 1 * (4) 

2t 7T t] 2 u t ^ u T 



— a-ir* 

u T J 



where r\ is the shear viscosity and (r w , X\) are second-order transport coefficients. 

In a DTT, the description of nonequilibrium hydrodynamic states requires the introduc- 
tion of a new tensor £ a7 w : 



evolution for ^ v reads 



rich is symmetric, traceless and vanishes in equilibrium [15|. The 



Al - U ^ Du f> (5) 



3r n r]u T M u 7 
and the shear tensor is calculated from the nonequilibrium tensor as follows 

W v = vr _ h^L^a^ _ Ia^M • (6) 

As independent variables we choose (p,u x ,u y ,U XX ,U xy ,U yy ) for the SOT and 
(p, u x , u y , £ xx , £ xy , E, yy ) for the DTT. The r component of the velocity follows from normaliza- 
tion, while the other nontrivial components of IP" (and of follow from the transversality 
and tracelessness conditions. We use vanishing initial transverse velocity and shear tensor, 
and set the initialization time r = 1 fm/c and the initial and freeze-out temperatures 
Tj = 333 MeV and Tf = 140 MeV, respectively. The initial energy density profile is calcu- 
lated using Glauber's model 35[ (see [8j for details). In all calculations we set the impact 
parameter 6 = 3 fm, and we use a 13 fm x 13 fm transverse plane and values for the second- 
order transport coefficients corresponding to a strongly-coupled M = 4 Super- Yang Mills 
(SYM) plasma : T n = 2(2 — ln2)r//(sT) and Ai = r]/(2nT), where s is the entropy 

density. In order to calculate the photon spectra, we take the critical temperature to be T c 
170 MeV which is consistent with recent Lattice QCD calculations (see 
therein) . 



and references 



The set of hydrodynamic equations must be closed with an equation of state. We employ 
two different EoS: the one computed by Laine and Schroder j^], which connects a high-order 
weak-coupling perturbative QCD calculation at high temperatures to a hadron resonance 
gas at low temperatures via an analytic crossover (as suggested by Lattice QCD calculations 



381]) . and an ideal one with p = p/3. In Fig. [TJwe reproduce the realistic EoS (which we will 
call LQCD EoS in what follows), showing the square of the velocity of sound c 2 s = dp/ dp 
as a function of temperature. We note that using the LQCD EoS and these values for r , 
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t w and Ai, data on Kaon elliptic flow, < px > anc 
reproduced both within the SOT and the DTT {2, 7. 



total multiplicity can be consistently 
SI- 




0.14 0.16 0.18 



FIG. 1: (Color online) Square of the velocity of sound (? a as a function of temperature corresponding 
to the EoS of Laine and Schroder [371] - This EoS connects a high-order weak-coupling perturbative 
QCD calculation at high temperatures to a hadron resonance gas at low temperatures via an 
analytic crossover. 



B. Photon production 



In this paper we consider thermal photon production due to the processes of Compton 
scattering, qq annihilation and bremsstrahlung in the QGP phase, and ttti — > pj, np — > 
and p — > 7T7T7 in the hadron phase. For a review on photon production in the context of 
heavy- ion collisions see 18 1. 

The Compton scattering and qq annihilation contribution is 



E 



dN 



d 4 xd 3 p 2tt 2 



(7) 
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where c ~ 0.23, a = 1/137 and [40] 

a s (T) 



6ir 



(8) 



(33 - 2N f ) ln(8T/T ( 

where Nf is the number of quark flavors. In the summation / stands for quark flavor while 
ef is the electric charge of the quark in units of the charge of the electron. The rate due to 
qq annihilation with an additional scattering in the thermalized medium is given by 

dN 



E- 



3tt 5 



aa s C£e 2 f )ETe- E / T (J T - J L ) 
f 



d A x(Pp 

where Jt ~ 1.11, Jl ~ 1.06. 

The rate of photon production due to bremsstrahlung reads 



(9) 



4l| 



E 



dN 
d A xd 3 p 



-aa. 



7P 



rp4 



~ E/T (Jt 



Jl 



3C(3) + 



7T 2 E 

~6T 



+ |J In 2 + m 3 (-e- E / T ) + 2^U 2 {-e- E ' T ) - ^ ln(l + e"^) 



(10) 



where ( is the zeta function and Li m = Yl^Li z " '/n m are polylogarithmic functions. 

In order to calculate the photon spectrum from the fireball when the LQCD EoS is used, 
one has to know also the photon production rate from the hot hadron gas since photons will 
also be emitted from this thermal phase following the QGP. We use the simple estimate for 
the photon production rate given by Steffen and Thoma [42], which reproduces the sum of 
production rates for the processes tttt — > p^y, np — > 7T7 and p — » txtv) calculated in 43] . The 
estimate for the production rate reads 

dN 



E- 



4gT 2.15 e -l/(1.35ET) 



0.77 _ E/T 



(11) 



d 4 xd 3 p 

The photon production rates just described are calculated with the equilibrium distri- 
bution function within Boltzmann's approximation, hence the factor e~ E ' T . However, it is 
well known that viscous effects lead to the modification of the thermal distribution func- 
tions which in turn modify observables such as hadronic elliptic flow [3l|. These 



modifications may also have significant observational effects on the photon spectra as has 
been shown recently 
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251 ]. It is therefore interesting to incorporate the viscous correction 



to the distribution function and asses its impact on final photon spectra. To this end, we 
calculate the photon production rates using the nonequilibrium distribution function coming 
from Grad's ansatz which is given by [l| 



/(.'^) = /o + /o(l-/o) 



2T2(p + p) 



(12) 



with fo the Fermi-Dirac function, and compare with the results obtained by using the equi- 

rat Grad's ansatz becomes unrealible at high 



librium distribution function. It is known t 
transverse momentum p T > 3— 4 GeV |g[ [25I 



44j and for this reason we will restrict to lower 



values of pt in what follows. 

The total thermal photon spectrum is obtained by integrating the sum of the production 
rates for each process given above over the space-time evolution of the fireball created in the 
collision. It is given by 

dN \ /• r* f Y dN 

dx / E - (13) 



d 2 pxdY J J J J_y d 4 xd 3 p 

where t\$, are the initial and final time one is interested in, Y is the rapidity of the nuclei, 
dx = (dx, dy) and pt is the transverse momentum. Note that the photon energy in the 
comoving frame is given by prcosh(y — Y'). 

III. RESULTS 

We now present our results on photon spectra, in all cases at midrapidity Y — 0. 

In Fig. [2] we show the total thermal photon spectra calculated with the SOT and the 
DTT for r)/s — 0.3, both for the ideal EoS and for the LQCD EoS, and including or not 
the nonequilibrium correction to the distribution function. We do not show the spectra 
calculated with the DTT for the ideal EoS because it is very similar to that obtained in the 
SOT. 

It is seen that the spectra strongly depend on the EoS used, being two orders of magnitude 
smaller in the case of the ideal EoS. This is a consequence of the fact that the hydrodynamic 
evolution of the plasma is slower when the realistic EoS is used and hence more thermal 



photons are produced (see [18l. 1231 125|). Besides, as already shown in |2l|, |25], the effect of 
the nonequilibrium correction to the distribution function is to harden the spectra. 

Going over to the comparison between the SOT and DTT, it is seen that the differences 
in the spectra arise entirely from the nonequilibrium corrections that become important 
only at > 1 GeV. The spectrum obtained from the DTT is significantly larger, and the 
difference in spectra increases with increasing p?. As we will show in what follows, this is 
due to the fact that the nonequilibrium corrections to the distribution function are larger 
in the DTT, and therefore the viscous correction to photon production is larger too (for a 
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FIG. 2: (Color online) Thermal photon spectra calculated with the SOT and the DTT for rj/s 
0.3, both for the ideal and the LQCD equations of state. 



discussion of this point in the context of hadronic elliptic flow see [8|). 

In Fig. |3] we compare the total photon spectra obtained with the LQCD EoS in the 
DTT and the SOT with rj/s = 0.3 and f]/s = 0.5. It is seen that the spectra is strongly 
dependent on the value of rj/s, in agreement with the results of 21(. Moreover, the difference 
in spectra corresponding to tj/s — 0.3 and rj/s = 0.5 is slightly larger in the DTT, showing 
that the latter formalism is more sensitive to changes in the value of the viscosity-to-entropy 
ratio. It is interesting to note that the spectra obtained in the SOT with rj/s = 0.5 is very 
similar to that obtained in the DTT with rj/s — 0.3, which clearly indicates that differences 
in hydrodynamic formalisms constitute a significant source of uncertainty in the extraction 
of rj/s from photon spectra data. A similar conclusion was found in [8| with respect to 
the charged hadron elliptic flow, although in that case differences between the SOT and 
the DTT have a smaller influence on the value of r]/s that can be extracted from data by 
matching to hydrodynamic results. 

In order to determine the influence of the nonequilibrium contribution on the different 
processes considered, in Figs. HHwe show this contribution integrated over the evolution 
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FIG. 3: (Color online) Thermal photon spectra calculated with the SOT and the DTT for Tj/s = 
0.3 and rj/s = 0.5. 

of the fireball in the SOT and the DTT for rj/s — 0.3, taking into account only the pro- 
cesses of gg-annihilation and Compton scattering, bremsstrahlung, or considering only the 
hadronic phase, respectively. It is seen that in all three cases the nonequilibrium correction 
is significantly larger in the DTT. The nonequilibrium correction to the spectra is larger in 
the hadronic phase at low pj>, as is the difference between the DTT and the SOT as it can 
be seen from Fig. |6j At larger values of px the nonequilibrium correction is larger in the 
QGP phase. The difference between the spectra calculated in the DTT and in the SOT is 
almost independent of Pt for the gg-anmhilation/Compton scattering and bremsstrahlung 
processes, while in the hadronic phase it decreases with increasing transverse momentum. 

The differences between the nonequilibrium contribution to photon spectra in the SOT 
and in the DTT can be understood from Fig. where we show the nonequilibrium distri- 
bution function of Eq. (fT2l) evaluated at (x, y) = and at r = 3 fm/c as a function of px 
for the DTT and the SOT with rj/s = 0.3 and Tj/s — 0.5. Several interesting features, some 
of which were already discussed when showing our results on spectra, can be seen from the 
figure. The difference between the DTT and the SOT becomes appreciable at pr > 0.5 GeV. 
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FIG. 4: (Color online) Nonequilibrium contribution to photon spectra calculated in the SOT and 
in the DTT taking into account only (/(/-annihilation and Compton scattering, for rj/s = 0.3. 

Beyond this value it is seen that the distribution function in the DTT is considerably larger 
than the one in the SOT. Moreover, the effect of increasing the value of rj/s is larger in the 
DTT, which reflects the fact that the viscous correction to f(x^,p^) is larger in the DTT. 

From these results we conclude that the nonequilibrium contribution to photon produc- 
tion strongly modifies the spectra, even at low transverse momenta, resulting in harder 



spectra (this is in line with the results of 2l|, |25|]). Moreover, we have found that the differ- 



ence between hydrodynamic formalisms leads to significant uncertainty in the extraction of 
7? Is from thermal photon spectra. 

nn 

It has been shown before [8|, |31[ that the nonequilibrium correction to the distribution 
function has also a strong influence on hadronic observables, especially on the elliptic flow. 
Therefore, the calculation of photon and hadronic observables using different hydrodynamic 
models may provide a way of determining or at least constrain, albeit phenomenologically, 
the form of the nonequilibrium correction to the distribution function. This may help 
improve the theoretical description of the freeze-out process, which is currently based on 
Grad's quadratic ansatz (see [3l|] for recent developments). 
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FIG. 5: (Color online) Nonequilibrium contribution to photon spectra calculated in the SOT and 
in the DTT taking into account only bremsstrahlung, for r//s = 0.3. 

IV. SUMMARY 



We have calculated the spectra of thermal photons produced in Au+Au collisions taking 
into account the processes of Compton scattering, gg-annihilation and bremsstrahlung in 
the QGP phase, and tttt — > fry, up — > tt7 and p — > tttt-j in the hadron phase. The calculation 
was done using two dissipative hydrodynamics formalisms to model the 2+1 evolution of the 
fireball, namely the second-order theory based on conformal invariance and a divergence-type 
theory. We have focused on determining the influence of the nonequilibrium correction to the 
distribution function on photon production in both formalisms. We have also compared the 
results obtained with an ideal and a realistic equation of state based on Lattice QCD that 
includes an analytical crossover between the quark-gluon plasma and a resonance hadron 
gas. 



In agreement with previous studies 



18 



23|-|25|) we have found that the use of a realistic 



equation of state enhances the photon spectra, implying that the QGP-hadron crossover can 
not be neglected in the investigation of thermal photons. Besides, the spectra calculated in 
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FIG. 6: (Color online) Nonequilibrium contribution to photon spectra calculated in the SOT and 
in the DTT taking into account only photon production in the hadronic phase, for rj/s = 0.3. 



the DTT is significantly enhanced with respect to the one calculated in the SOT. The dif- 
ference is entirely due to differences in the nonequilibrium corrections to photon production 
in both formalisms. 

Our results point to the conclusion that differences in dissipative hydrodynamic for- 
malisms are a significant source of uncertainty in the precise determination of rj/s from 
data, particularly from those observables, such as the one studied here, that depend on the 
nonequilibrium distribution function. Comparing our results with those obtained for the 
charged hadron elliptic flow and its matching to RHIC data jsl, we find that the differ- 
ence in hydrodynamic formalisms used to calculate thermal photon spectra leads to a larger 
uncertainty in the possible extraction of rj/s from measured photon spectra. 

It is interesting to note that the dependence of hadron or photon observables on the 
nonequilibrium correction to the distribution function, which leads to uncertainty in the 
extraction of rj/s from data, may in turn serve to constrain the form of the nonequilibrium 
distribution function. In this way one could gain some insight on the relation between 
kinetic theory and hydrodynamic models as well as on possible generalizations of Grad's 
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FIG. 7: (Color online) Nonequilibrium distribution function evaluated at the center of the fireball 
and at r = 3 fm/c as a function of px calculated in the DTT and the SOT with rj/s = 0.3 and rj/s 
= 0.5. 

quadratic ansatz, and help improve the theoretical description of the freeze-out process in 



heavy-ion collisions [31|. For this to be feasible, the theoretical uncertainties present in the 
hydrodynamic description of heavy-ion collisions, for example those regarding the initial 
conditions, must be under control. 

In this paper we have limited ourselves to the case of vanishing bulk viscosity. As men- 
tioned in the Introduction, Lattice QCD results indicate that the bulk viscosity becomes 
peaked at the critical temperature. It would be interesting to include bulk viscosity in the 
DTT to determine its influence on thermal photon spectra in this hydrodynamic model. 
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